Non-gaussianity in axion N-flation models 
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We study perturbations in tlie multi-field axion N-flation model, taking account of the full cosine 
potential. We find significant differences with previous analyses which made a quadratic approxima- 
tion to the potential. The tensor-to-scalar ratio and the scalar spectral index move to lower values, 
which nevertheless provide an acceptable fit to observation. Most significantly, we find that the 
bispectrum non-gaussianity parameter /nl rnay be large, typically of order 10 for moderate values 
of the axion decay constant, increasing to of order 100 for decay constants slightly smaller than the 
Planck scale. Such a non-gaussian fraction is detectable. We argue that this property is generic in 
multi-field models of hilltop inflation. 

PACS numbers: 98.80.Cq 



Much focus has been placed lately on the discov- 
ery potential of cosmic non-gaussianity in the statistics 
of primordial perturbations. The Wilkinson Microwave 
Anisotropy Probe (WMAP) has already set interesting 
limits |l|. The Planck satellite, now taking data, will 
improve these significantly, reaching a sensitivity to the 
non-gaussianity parameter /nl of around five. Discovery 
of non-gaussianity would open a new arena of cosmo- 
logical observations particularly suited to probing early 
Universe physics. 

Present ideas in fundamental physics suggest there 
may be many scalar fields which can influence the early 
Universe, including inflation. N-flation (2| uses many 
string axions to provide a realization of the 'assisted in- 
flation' phenomenon [3i]i in which a collection of scalar 
fields cooperatively support inflation even if their poten- 
tials are individually too steep. The phenomenology of 
such models (see also Ref. [4|) links fundamental physics 
and upcoming cosmological observations. 

Previous N-flation studies have assumed that all rele- 
vant fields are close to their minima and can be described 
by quadratic potentials. For axions the full potential is 
trigonometric and we find the quadratic approximation 
is unreliable. Even for identical potentials, the condition 
for stable co-evolution of the fields is violated near the 
hilltop y. Therefore fields in this region evolve on di- 
vergent trajectories. Accounting for this divergence by 
retaining the full potential leads to two very significant 
changes. The predicted scalar spectral index and tensor- 
to-scalar ratio, r, are reduced. This remains compatible 
with existing observations but may leave r undetectable. 
More importantly, /nl is predicted to be large, and very 
plausibly within the range of future probes. 

This unexpectedly large non-gaussianity is a genuine 
multi-field phenomenon. It is a consequence of the di- 
verging trajectories near the hilltop, implied by a nega- 
tive ?7-parameter of order unity or larger. In single-field 
models, potentials of this form lead to a density pertur- 
bation with a spectral index, n, in conflict with obser- 



vation. The assisted inflation mechanism reduces 1 — n 
to an acceptable value, but leaves /nl dominated by the 
contribution of the field closest to the peak. 

THE MODEL 

The axion N-flation model is based on a set of A^f un- 
coupled fields, labelled 0i, each with a potential [3] 



V^ = Af (1-cosa,) 



(1) 



where a^ — liK^ij ji and /j is the i axion decay con- 
stant. In a more general model couplings may exist be- 
tween the fields, but we will not consider these. The 
mass of each field in vacuum satisfies mi — 2itK^J fi, and 
the angular field variables ai lie in the range (— tt, +Tr]. 
Without loss of generality we will set initial conditions 
with all Ui positive. If only a single field is present this 
model is known as natural inflation [6|. 

Calculation of the observables n, r and /nl makes use 
of the 5N formula [7| , which considers how the total num- 
ber of e-folds of expansion A^ is modified by field pertur- 
bations. We define slow-roll parameters for each field as 



Ml 



(2) 



where Mp = (SttG)^^/^ is the reduced Planck mass, a 
prime denotes the derivative of a function with respect 
to its argument, and no summation over i is implied. 
The global slow-roll parameter e = —H/H^ can be writ- 
ten as a weighted sum e ~ X)i(^i/^)^^i! i^ which each 
field contributes according to its share of the total energy 
density. We must have e < 1 during inflation. 

We work in the horizon-crossing approximation, in 
which the dominant contribution to each observable is 
assumed to arise from fluctuations present only a few 
e-folds after horizon exit of the wavenumber under dis- 
cussion. After smoothing the universe on a superhori- 
zon scale somewhat smaller than any scale of interest. 



the horizon-crossing approximation becomes valid when- 
ever the ensemble of trajectories followed by smoothed 
patches of the universe approaches an attractor. We sup- 
pose that inflation exits gracefully, with each field settling 
into the minimum of its potential. The horizon-crossing 
formulas will then be a reasonable approximation. Us- 
ing this method, and conventional definitions for each 
observable parameter [81] , we find 
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where N ,i and A^^y are respectively the first and second 
derivatives of TV with respect to the fields, and * indicates 
evaluation at horizon crossing (determined by Eq. ^ be- 
low). In writing Eq. (J6|) any intrinsic non-gaussianity 
among the field perturbations at horizon crossing has 
been neglected, a good approximation provided /nl > 1 
[3, [ifl] • Our sign convention for /nl matches WMAP [1| , 
and the non-gaussianity is predicted to be of local type. 
The observed amplitude of perturbations is obtained by 
adjusting the A^ to give an appropriate value of 7J, . 

Under a quadratic approximation to each potential, it 
can be shown that Eqs. ^ and Kl recover their single- 
field values of order ~ 1/A^* [idlulj rnaking /nl unde- 
tectably small. The spectral index can be shown to be 
less than its single field value 1 — 2/iV, [12| with equality 
only in the equal-mass case. Its value for a given choice 
of parameters must be computed numerically [13| . How- 
ever, we will see that these results all change whenever 
our initial conditions populate the hilltop region. 



N-FLATION PERTURBATIONS 

Eqs. ©-(O apply for any choice of A^ and ji. We re- 
strict attention to the case where all fields have the same 
potential, which already captures the interesting phe- 
nomenology. A broader investigation will be published 
elsewhere. The scale A = A^ is fixed from the observed 
amplitude of Vq, leaving f = fi and Nf as adjustable 
parameters. The initial conditions are drawn randomly 
from a uniform distribution of angles a^ , with several re- 
alizations to explore the probabilistic spread. From these 
two parameters we predict the observables n, r and /nl- 

There are two constraints. First, we require sufficient 
e-foldings. For a given set of initial angles a^, and ig- 
noring a small correction from the location of the end of 
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where in the second equality we have replaced A'tot with 
its expectation value by averaging over a^. Eq. ([7]) is 
replicated to high accuracy in numerical simulations. For 
a given / it determines the minimum number of fields re- 
quired for sufficient inflation, typically several hundred or 
more. There is no similar constraint from the spectral in- 
dex. When A'tot ~ A^*, the a* are uniformly distributed 
and (n. — 1} — —5 In2/A',, independent of / and A^f. This 
tilt is observationally acceptable. For larger A'f the spec- 
tral index approximately satisfles Eq. (|5]) below. 

Second, a key motivation of the N-flation model was to 
obviate the requirement for superplanckian field values, 
which are invoked in many single-field models. If one 
literally imposes \(^\ < Mp this requires fi < 2Mp for 
each i. However, it would be reasonable to regard this 
condition as approximate and not mandatory. 

The ei approach zero for fields close to the hilltop, so 
each summation in Eqs. ((21)^® is dominated by those 
fields with the smallest e^. Suppose some number N of 
fields have roughly comparable e^, of order e. The ob- 
servable parameters have different scalings with A'^. The 
spectrum, V( , scales like N copies of a single-field model 
with slow-roll parameter e, whereas r is reduced by a 
factor N compared to its value in the same single-field 
model. The spectral index can be written exactly (within 
slow-roll) in terms of a single sum coming from 77* , 

n-l^-2e.-8n'(^^^ /^(1-cosa*), (») 

and is independent of N. It becomes close to — 2e, when 
the denominator is of order 10'^. This is the standard 
assisted-inflation mechanism. Most importantly, /nl has 
the approximate behaviour 
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(9) 



which is independent of e if the dominant fields are suffi- 
ciently close to the hilltop. N-flation has lifted the single- 
field consistency condition /nl ~ ^(5/12)(n— 1) [9|,|l0|, 
which prevents single-field models generating large non- 
gaussianity without violating observational bounds on n. 
Where the summations in Eqs. ©-(11]) are dominated 
by a single field, this formula shows that /nl can become 
rather large, scaling as (Afp//)^. For / = Mp, we find 
/nl ^ 16.4; a non-gaussian fraction of this magnitude 
should be visible to the Planck satellite. It is even possi- 
ble to achieve /nl ~ 100 for / ~ 0.4Mp, though then A^f 
must be very large to gain sufficient e-foldings. If /nl ^ 
50 it may be more profitable for Planck to search for non- 
linearity in the trispectrum 1141], for which estimates in 
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FIG. 1: Predictions in the n-r plane, averaged over realiza- 
tions, for various values of / between 0.4Mp and 2Mp and of 
N{ between 464 and 10,000, all giving sufficient inflation. The 
black (left) cluster of points takes A^, = 50 and the red (right) 
cluster A^* = 60. The quadratic expansion predicts r = 8/N^, , 
far off the top of this plot. The region right of the line is 
within the WMAPT+BAO+ii'o 95% confidence contour [I[. 



FIG. 2: Predicted non-gaussianity, |/nl, for / — Mp and 
Af* = 50. The error bars are on the mean over realizations 
(not the standard deviation). Here the maximum achievable 
value of |/nl is 2-k'^ ~ 20, almost saturated in some real- 
izations. The significant spread is due to initial condition 
randomness with typical mean values being around half the 
maximum achievable value, and no discernable trend with A^f . 



the quadratic approximation were given in Ref. jl5[ . We 
defer a full analysis of the trispectrum to future work but 
note that the trispectrum equivalents of Eq. ^ are, in 
conventional notatioii [H, tnl = (47r'*/iV2)(Af4/j4) ^^^ 
(54/25)5NL = (87rViV2)(M4//4). 

The expectations described above are borne out in nu- 
merical calculations. In Fig. [1] we show model predictions 
in the n~r plane, averaged over several realizations of the 
initial conditions. We see n and r are only weakly depen- 
dent on the model parameters (though there is significant 
dispersion amongst realizations, not shown here), with 
the choice of N^ being the principal determinant of n. 
In Fig. [2] we plot /nl as a function of N{ for / = Mp, 
with ten realizations at each A^f. This clearly shows the 
expected maximum, which is nearly saturated in cases 
where a single field dominates the summations. In cases 
where several fields contribute significantly to the sums in 
Eqs. (I3|)-(l6|), the non-gaussian fraction is reduced. Fig. [3] 
shows the mean predicted non-gaussianity, averaged over 
realizations, as a function of /. 

Eqs. ^ and ^ clarify the origin of large /nl in this 
model. The cooperative effect of the N-fiation mecha- 
nism does not enhance the non-gaussian signal. Indeed, 
/nl is suppressed by the central limit theorem where 
N :^ 1 fluctuations contribute equally to the curvature 
perturbation. Nor does the large effect arise from a sin- 
gularity in the e-folding history, A^, as a function of its 
initial angles ai. Although Eq. ([7]) is singular in the 
limit ai -^ tt, its Taylor expansion is trustworthy un- 
less \ai — 7r| < (7'^r)^/^(Mp//i). The observed magni- 
tude of V^ requires \ai — 7r| > r^/^(/i/Mp) for each field, 
so a breakdown of the Taylor expansion cannot become 



relevant unless at least one fi is a few orders of magni- 
tude less than the Planck scale, of order (/^/Mp)'* < P(^. 
These constraints additionally imply that we do not tres- 
pass on any region of field space where quantum diffusion 
competes with classical motion. 

Instead, the large /nl derives from a generic dispersive 
effect present in any hilltop potential. Measuring the 
displacement of 4>i from the hilltop by Si , each potential 
can be approximated in its vicinity by Vi ~ 2A^(1 + 
r]iSf/2Mp), where 77.^ < satisfies 
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These potentials are tachyonic. Fields close to the hill- 
top remain almost stationary, while fields further away 
are ejected downhill. This process typically leaves a few 
fields on top of the hill, which have small e^ and dominate 
the sums in Eqs. ([3])-([ni)- It seems clear this behaviour 
is generic for any N-flation model constructed using hill- 
top potentials. The few fields remaining in the vicinity 
of the hilltop each generate contributions to the curva- 
ture perturbation with third moment (6/5)/nl ~ —11* |9|. 
Accounting for suppression arising from the central limit 
theorem, we recover the approximate expression (|9]). For 
a general hilltop potential, well-rehearsed arguments lead 
us to expect |??| ~ 1 and therefore /nl "^1- In a single- 
field model this is the 'ry problem'. In an N-flation model, 
it is a generic expectation of enhanced non-gaussianity. 
Even larger yields are possible in some models, includ- 
ing our case, if it is possible to achieve I??! 3> 1 while 
preserving technical naturalness. 
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FIG. 3: The predicted non-gaussianity as a function of /, for a 
range of choices of A'f . Each point shown is the average of five 
or more realizations for an f~N{ pair. We see a strong trend 
with /, well represented by Eq. ((9} with TV ~ 2. The different 
A^f are scattered by randomness in the initial conditions rather 
than an identifiable trend. 



CONCLUSIONS 

We have described a new mechanism for generating 
observably-large cosmic non-gaussianity, based on the 
strongly-dispersive dynamics of fields in a hilltop region. 
In a multi-field context such as the axion N-flation model, 
assisted inflation can yield a viable spectral index with- 
out a major dilution of the non-gaussianity. As compared 
to the quadratic potential approximation to N-flation, we 
found a substantial decrease in r, a modest increase in 
1 — n, and a substantial increase in /nl- These changes 
will happen whenever initial conditions have a significant 
probability of populating the hilltop region, such as the 
uniform (in field angle) initial conditions we chose. 

Searches have previously been made for models which 
achieve | /nl | ^ 1 while preserving slow-roll during in- 
flation [17|. The N-flation model is of this type, but 
offers several advantages. The non-gaussian fraction is 
naturally bounded above, so that /nl cannot become ar- 
bitrarily large. Therefore our predictions do not depend 
on a sudden exit from inflation, e.g. triggered by a hybrid 
transition, to prevent /nl from growing to an unaccept- 
able value. Equally important, our large signal does not 
derive from a singularity of the e-folding history iV, as a 
function of its initial conditions. These means we can rely 
on a perturbative expansion. We can simultaneously sat- 
isfy observational constraints on the spectral index and 
tensor fraction. Moreover, this result seems generic. In- 
flation is self-replicating on top of the hill, sometimes de- 
scribed as 'topological inflation' [18|. Coupled with the 
dispersion of trajectories originating from the vicinity of 
the hilltop, this implies large non-gaussianity may not be 
uncommon over a landscape of scalar field vacua. 
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